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Abstract 

Bacterial colonies, granular fingers, and surfaces generated by point invasion all involve circular 
patterns. As the pattern grows, the size of interface increases. And the Kardar-Parisi-Zhang and 
Edwards- Wilkinson equations break down because they all fail to conserve the fluid volume for a 
moving boundary. If we plot the interface against the polar angle, the wriggling of the interface 
wrinkles progressively with time. It shows convincingly that the noise spectrum behind the growth 
of a circular pattern is drastically different from the white noise commonly adopted for a planar 
interface. This ruins the universal power law found in the latter system. A new model was 
proposed for our revisit of the Hele-Shaw experiment using a circular invasion. Analytic solutions 
are obtained in excellent agreement with experimental results. 

PACS numbers: 68.35.Ct, 68.35.Fx, 68.35.Rh, 46.65. +g 
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The properties related to surface growth have attracted much attention. In the past 
wo decades, scientists reported scaling behavior in many experiments, like paper wetting 
, burn front [3|, and fluid invasion into disordered media etc.. Surprisingly the 

surface roughness of systems with different growth mechanisms turns out to share the same 
scaling law. In 1989, Rubio, Edwards, Dougherty, and Gollub found the scaling behavior to 
obey a power law with exponent a, w ~ L a J4]. Similar to the critical phenomena at phase 
transitions, this exponent is used to distinguish different universality classes and character- 
ize their scaling behavior Among the various models for the growth mechanism, the 
Kardar-Parisi- Zhang (KPZ) model is perhaps the most famous 8J. It includes a nonlinear 
term to account for the askew growth in the diffusion equation of the Edwards- Wilkinson 
(EW) model jjj. Several papers have tried to introduce more accurate terms to describe 
the depinning dynamics of a driven interface in an isotropic disordered media, such as the 
quenched EW jsj], quenched KPZ {§], and quenched Herring- Mullins for a tensionless inter- 
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111 ]. However, discrepancies are always there when compared to the experiments 



(See Tab. I). So we decide to revisit the classic Hele-Shaw experiment. 

We choose to inject the fluid from a circular inlet in contrast to the standard planar 
one. Similar experiments do exist in the literature, but they often involve complications 
like the use of non-Newtonian fluid or interface instabilities by prefilling the media with 



another immiscible fluid [12j, |13j. Naively, we might expect the surface generated by a 
circular injection to eventually approach that of a planar one as the area of the invasion 
liquid becomes large. We emphasize in this Letter that this presumption is wrong mainly 
because the circumference of the circularly injected fluid grows with time. This causes the 
constant introduction of new local structures as the surface expands, while old ones being 
pushed apart but, nevertheless, reserved (Fig.l). 

Since the area of fluid is J R 2 d8/2 for our circular invasion as opposed to J hdx for the 
planar case, R 2 is a more natural quantity to define the interface width than R. We first 
expand it in Fourier series: 

00 

R 2 (0, t) = R 2 (t) + \ a i (*) cos 19 + b i (*) sin W \ W 
1=1 

where R$ is the mean radius of the whole interface. Then substitute it into the definition of 
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interface width, 
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where AO denotes the range of the angular average, and iV = 2tt/A6 is the number of 
divisions of the interface. It is easy to show that the first term in the square root equals 
-Rq (t). The second term can be simplified if the roughness is small enough (ai,bi <C R^) 
to warrant the Taylor expansion. In contrast to the random term in the KPZ model, 
we include the randomness by ensemble-averaging over the Fourier coefficients of Eq.(l) 
in the final expression for the interface width. Physically, this randomness comes from 
the distribution of spacing between neighboring beads along the advancing front. Overall, 
narrower spacings enjoy a stronger capillary force and thus a higher growth rate. Note that, 
in contrast to the white noise for a parallel injection which is uncorrelated in both sites and 
heights, the increasing size of the interface introduces chronologically high Fourier modes 
into the surface (Fig.l), which forces us to adopt a decreasing /—dependence for the noise, 
(aidi') = (bibi>) = (?5ij' and (ai) = (bi) = where the angular bracket denotes taking the 
statistical average. 

This allows us to drop most of the cross terms in the double summation over I and keep 
only the relevant terms. The final expression for the width becomes 
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Schematic plot of this result can be found in Fig. 2. Note that when AO is small, this width 
always scales linearly with AO and gives a = 1. In fact, this conclusion does not rely on the 
nature of noise spectrum. Without it, the equivalent expression of Eq.(3) becomes messier, 
but the lowest order term is still linear in AO when Taylor expanded. Since the next order 
term has a negative coefficient, it becomes natural when we try to fit AO — c (AO) 2 by (A0) a 
that an exponent both less than 1 [4J, [5|, laj and time-dependent is obtained. So the eventual 
magnitude of this exponent becomes arbitrary and varies with the range of AO we adopt for 
the power law fitting. This therefore casts doubts on the scaling arguments people conjure 
up to derive a and conclude that it plays the same universal role as the exponents in the 



critical phenomena 
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In order to study the dynamical behavior of Eq.(3), the actual dependence of £ on / 
is required. Note that our system size is growing in time. This is a moving boundary 
problem, and both the EW and KPZ equations fail to conserve the area. Instead of further 
complicating the existing models, we adopt the polar coordinates and modify the EW model 
to guarantee this conservation. 

^ = "-^ + t(t)[R 2 (t,e)-RUt)] + ^ (3) 

where V a denotes the speed of fluid injection. Two phenomenological coefficients are in- 
troduced. First, the relaxation coefficient v is the same as in the EW and KPZ models, 
which derives from the surface tension and favors a smooth interface. The second term on 
the RHS is new. We ascribe it to the inertia effect which tends to magnify the existing 
crests by allowing them to advance faster in the next moment than the troughs. Therefore, 
we expect its coefficient £ to be proportional to the average kinetic energy of the front, 
cr(dRo/dt) 2 ~ crVa/t. It is interesting that Eq.(4), under the change of variable from 9 
to x = R • 9, will generate a nonlinear term (dR/dx) 2 characteristic of the KPZ model. 
Dimensional analysis allows us to pinpoint the dependence of these two phenomenological 
coefficients on the microscopic parameters as: £ (t) — £ jt — KaV a /r]r 2 t and v = Tr/rj where 
a, ft, r], r, T are the fluid density, permeability, viscosity, bead size, and dynamic surface ten- 
sion. The reason why the contact angle did not come in is that we have limited the spacing 
between the two acrylic plates to be just a few beads to minimize its role. Note that the 
source term in Eq.(4) vanishes when we plug Eq.(l) into Eq.(4). Applying the orthonormal 
relation, we end up with a simple linear PDE for each of a/ and hi. They can be easily solved 
to determine the evolution of the average noise spectrum Q as, 

Cf(t) = Cf(t )e-^ 2 ^\t/t )^ (4) 

Substituting the new noise form into Eq.(3), we obtain 



w = m. 



^g(i-w )e " W(, " ,0,c ' 2(i ° ) (5) 

ZKo to 

where N' is the fluid-covered area divided by glass bead cross-section and (i(to) denotes the 
initial noise spectrum at to- Since there is no explicit random term in our PDE, this initial 



pattern cannot be taken as a perfect circle, otherwise the subsequent evolution will not be 
able to produce any roughness. But, we also need to make sure later that the final theory 
does not depend on the specific choice of this initial time t and its corresponding pattern 
Ci(to)- We considered several possible situations and generated the spectrum of the noise 
by analyzing the interface. The final form we adopt to describe the decay of noise with / 
is Cz(io) = + / + s'l 2 ) where the value of s is found to fall between 0.1 and 1, and s' is 
around 0.001. This corresponds to an initial shape in Fig.2. 

Our experimental setup includes a thin cell (0.15 cm deep) formed by a horizontally 
positioned acrylic dish (radius=17 cm) and an acrylic top plate, a fluid injection system, 
and a data recorder. The syringe injects the fluid vertically into the disorder media from 
the center of the dish with speeds ranging from 2.19 to 2.50 cm 2 /s. Glassbeads of diameter 
500 fim are well tamped between the acrylic plates to make sure they will not move during 
the invasion. We use a high resolution (1 million pixels) CCD camera to monitor and record 
the interface, which is set vertically above the dish center. The camera monitors an area 
of 24 x 24 cm 2 plus the edge space, and the size of each pixel is about 450 fim 2 . After the 
picture has been taken, a personal computer digitizes the video image and calculates the 
interface width. 

Experimentally the weighting of the noise spectrum is found to decrease with the Fourier 
angular momentum. This is evidenced by Fig.l, which shows the interface wiggles progres- 
sively more severe with time. The high modes come in as more local structures sprout up 
at the newly gained surface. Lower modes, which correspond to more global structures, 
emerge earlier and thus have a longer history to evolve. According to the instability term in 
Eq.(4), fluctuations are amplified via feedbacks on themselves. Therefore, older structures 
tend to have larger amplitudes. The scale dependence of the width is shown in Fig. [5j The 
linear behavior of the width with a slope exactly equal to one in the asymptotic regime of 
small scales is consistent with we expected. Same features can be found in the bulk of non 
super- rough interface and many recent experiments 0, Q, 19] . 



We use the parameters in Tab. II to estimate £' and v at about 1 and 0.036. Experimental 
values put them at about 0.8 ~ 1.1 (inset of Fig.4) and 10~ 2 ~ 10 _1 respectively. Other 
features are also verified experimentally, which include: (1) The instability coefficient £' 
is linearly proportional to the invasion velocity (Fig.4), (2) no dependence on V a for the 
relaxation coefficient v. Moreover, Fig. 5 finds A6 dependence for j3, as has been reported 
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for a parallel injection^]. By use of Eq.(6), we can obtain the initial noise spectrum Q(to) 
by fitting the scale dependence of the interface width. Nicely, we deduce the same form 
at different times. Note that it is necessary to incorporate the possible symmetry of the 
initial configuration before being able to determine it directly from the spectrum in Fig.2. 
We duly checked that the interface width was not sensitive to the exact choice of this initial 
time to and shape of interface as long as they belong to the early stage when the effect 
of randomness is still small. As illustrated above, both a and (3 are smooth and varying 
functions of the scale. On the other hand, £' and v are good parameters to describe the 
scaling and dynamic behavior of our circular interface. 

Note that with our previous choice of a constant injection speed the advancing interface 
is decelerating. In order to check that our conclusions are not its artifacts, we repeated by 
accelerating the injection flux to attain a constant interface velocity. Now, £ = anV^/rjr 2 
becomes time- independent, V a is replaced by V^t, but parameter v is the same. The ex- 
perimental results in Fig.6-(a),(b) turn out to be qualitatively the same as Fig.3. We find 
that £ is still proportional to V' a with a slope that falls rightly within our estimated range 
an/fir 2 = 10 3 ~ 10 4 (Fig.6-(c)) and v remains independent of V' a with the right magnitude 
(Fig.6-(d)). We believe that if the value of k can be more accurately measured, the correct- 
ness of an/rir 2 shall be more solidly verified. We can now safely say that our dimensional 
analysis in the constant flow rate remains valid for the accelerated flow rate. 

In conclusion, we pointed out that the roughness of a circular interface is intrinsically 
different from the standard planar injection. Major difference lies in the fact that the former 
interface grows in size. This results in a noise spectrum whose amplitudes decrease with 
momentum, in contrast to the white noise for a planar interface. One major consequence 
is the lack of power law in the scaling and dynamic behavior, which we believe should be 
common to all circular patterns 20[]. Interestingly, we also find that the time exponent varies 
with the range of scale, as has been reported [6|] for a planar interface but not captured by 
existing models. All these features can be described by a linearized model we modified 
from the EW equation. Two phenomenological coefficients are put forward to represent 
the competing trend of instability and smoothness. Their dependence on the microscopic 
parameters are determined from dimensional arguments, and estimated values conform with 
the experiments. 

We benefit from discussions with Profs. H. H. Lin and C. K. Chan, and thank Prof. W. 
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TABLE I: Comparison of the scaling exponent in the literature using a planar injection 



Experiment a Model a 

Rubio[4] 0.73 KPZ[8] 0.5 

He[5j 0.65~0.91 RFIM[14] 1 

Soriano[6j 0.94 TOB[15j/TBO[16] 0.7/0.6-0.9 



TABLE II: The values of the relevant parameters in the experiment at room temperature. 



Parameter 


Value (unit) 


water-air surface tension (T) 


0.0727 N/m 


dynamic viscosity (77) 


1.002 x 10" 3 Pa-s 


permeability (re) 


~ 10~ 9 m 2 



160 


Saturated:17s 


i«o 

'a. 120 




» 100 


Transient:!}. 7s 


TO 80 

or 

60 


Early time:2.5s 


-3-2-10123 
9 (rad) 



FIG. 1: Time-lapse interface profiles for invasion speed of 3.31 cm 2 /s. More local structures emerge 
with time. 




FIG. 2: Schematic of the initial configuration of the invading fluid reconstructed from 0(irj) 
1 (0.2 + I + 0.005Z 2 ), obtained from the solid-line fit of the actual roughness spectrum Ci(to) with 
an invasion velocity of 5.1cm 2 /s at transient times in the inset. 



8 




FIG. 3: Logarithmic plot of the time dependence of m defined in Eq.(9) for three invasion speeds, 
V a \ = 2.50, V a 2 = 2.40, and V a 3 = 2.19cm 2 /s. The solid line in the inset indeed passes through 
the origin when extrapolated. 



12h 0.6f 




FIG. 4: Interface width versus time for invasion velocity 3.15 cm 2 /s under three different AO. The 
forms of these three solid lines are determined by the model. The inset illustrate (3 against AO for 
invasion velocity 3.15 cm 2 /s. 
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0.4 o 5 6 Saturated:25s 
3(rad) 




Transient:14s 



Early time:2.5s 



2 3 4 5 
AG (rad) 



FIG. 5: Interface width versus AO at three different stages labelled by early time, transient, and 
saturation. The solid line comes from our formula. Invasion velocity is 5.1cm 2 /s. The inset 
demonstrates that the actual value of a varies with the range of AO we use to do the power law 
fitting. 
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FIG. 6: (a) Scale dependence of a constant advancing interface with V' a = 4.01pixel 2 /s 2 . Solid line 
is the theoretical prediction. Because the flow rate is slow initially, more time is required for the 
designated early stage than in Fig.4. In (b), a is shown to vary with the scale, (c) £ is plotted 
against V' a . The slope of the fitting line, 6840, falls within our predicted range, (d) Experimental 
values of v are independent of V^, consistent with our theory. 
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Circular invasion of fluid into a quenched disordered media 
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Department of Physics, National Tsing Hua University, Hsinchu 30043, Taiwan, Republic of China 
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Bacterial colonies, granular fingers, and surfaces generated by point invasion all involve circular 
patterns. As the pattern grows, the size of interface increases. And the Kardar-Parisi-Zhang and 
Edwards- Wilkinson equations break down because they all fail to conserve the fluid volume for a 
moving boundary. If we plot the interface against the polar angle, the wriggling of the interface 
wrinkles progressively with time. It shows convincingly that the noise spectrum behind the growth 
of a circular pattern is drastically different from the white noise commonly adopted for a planar 
interface. This ruins the universal power law found in the latter system. A new model was proposed 
for our revisit of the Hele-Shaw experiment using a circular invasion. Analytic solutions are obtained 
in excellent agreement with experimental results. 

PACS numbers: 68.35.Ct, 68.35.Fx, 68.35.Rh, 46.65.+g 



The properties related to surface growth have attracted 
much attention. In the past two decades, scientists re- 
ported scaling behavior in many experiments, like pa- 
per wetting j burn front Q , and fluid invasion into 
disordered media 0, HI, [|[ etc.. Surprisingly the surface 
roughness of systems with different growth mechanisms 
turns out to share the same scaling law. In 1989, Ru- 
bio, Edwards, Dougherty, and Gollub found the scaling 
behavior to obey a power law with exponent a, w ~ L a 
Similar to the critical phenomena at phase transi- 
tions, this exponent is used to distinguish different uni- 
versality classes and characterize their scaling behavior 
Q. Among the various models for the growth mecha- 
nism, the Kardar-Parisi-Zhang (KPZ) model is perhaps 
the most famous It includes a nonlinear term to 
account for the askew growth in the diffusion equation 
of the Edwards- Wilkinson (EW) model @ . Several pa- 
pers have tried to introduce more accurate terms to de- 
scribe the depinning dynamics of a driven interface in 
an isotropic disordered media, such as the quenched EW 
Q, quenched KPZ Q, and quenched Herring-Mullins for 
a tensionless interface [T3, El- However, discrepancies 
are always there when compared to the experiments (See 
Tab. I). So we decide to revisit the classic Hele-Shaw ex- 
periment. 

We choose to inject the fluid from a circular inlet in 
contrast to the standard planar one. Similar experiments 
do exist in the literature, but they often involve compli- 
cations like the use of non-Newtonian fluid or interface 
instabilities by prefilling the media with another immis- 
cible fluid [l2|, UM ■ Naively, we might expect the surface 



TABLE I: Comparison of the scaling exponent in the litera- 
ture using a planar injection 



Experiment 


Q 


Model 


a 


Rubio[4j 


0.73 


KPZ[8j 


0.5 


He[5] 


0.65~0.91 


RFIM[14] 


1 


Soriano[6] 


0.94 


TOB[15]/TBO[M] 


0.7/0.6^0.9 
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FIG. 1: Time-lapse interface profiles for invasion speed of 
3.31cm 2 /s. More local structures emerge with time. 



generated by a circular injection to eventually approach 
that of a planar one as the area of the invasion liquid be- 
comes large. We emphasize in this Letter that this pre- 
sumption is wrong mainly because the circumference of 
the circularly injected fluid grows with time. This causes 
the constant introduction of new local structures as the 
surface expands, while old ones being pushed apart but, 
nevertheless, reserved (Fig.l). 

Since the area of fluid is J R 2 d0/2 for our circular in- 
vasion as opposed to J hdx for the planar case, R 2 is a 
more natural quantity to define the interface width than 
R. We first expand it in Fourier series: 



R 2 (6,t)=R 2 Q (t) 



OO 

E 

i=i 



ai (i) cos W + bi (t) sin 10] (1) 



where Rq is the mean radius of the whole interface. Then 
substitute it into the definition of interface width, 
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where A8 denotes the range of the angular average, and 
N = In j AO is the number of divisions of the inter- 
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Invasion center 



The cavity between 
porous media 



FIG. 2: Schematic of the initial configuration of the invading 
fluid reconstructed from Ci(io) = l/(0.2+/+0.005/ 2 ), obtained 
from the solid-line fit of the actual roughness spectrum C((io) 
with an invasion velocity of 5.1 cm 2 /s at transient times in 
the inset. 



face. It is easy to show that the first term in the square 
root equals R^{t). The second term can be simplified 
if the roughness is small enough (ai, bi <C Rq) to war- 
rant the Taylor expansion. In contrast to the random 
term in the KPZ model, we include the randomness by 
ensemble- averaging over the Fourier coefficients of Eq.(l) 
in the final expression for the interface width. Phys- 
ically, this randomness comes from the distribution of 
spacing between neighboring beads along the advancing 
front. Overall, narrower spacings enjoy a stronger cap- 
illary force and thus a higher growth rate. Note that, 
in contrast to the white noise for a parallel injection 
which is uncorrelated in both sites and heights, the in- 
creasing size of the interface introduces chronologically 
high Fourier modes into the surface (Fig.l), which forces 
us to adopt a decreasing I— dependence for the noise, 
{aia v ) = (bibi>) = (?5i,i< and (a t ) = (bi) = where the 
angular bracket denotes taking the statistical average. 

This allows us to drop most of the cross terms in 
the double summation over I and keep only the relevant 
terms. The final expression for the width becomes 
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Schematic plot of this result can be found in Fig. 2. Note 
that when A9 is small, this width always scales linearly 
with AO and gives a = 1. In fact, this conclusion does 
not rely on the nature of noise spectrum. Without it, 
the equivalent expression of Eq.(3) becomes messier, but 
the lowest order term is still linear in A8 when Tay- 
lor expanded. Since the next order term has a nega- 
tive coefficient, it becomes natural when we try to fit 
AO — c(A9) 2 by (A9) a that an exponent both less than 
1 @, H, Q and time-dependent is obtained. So the even- 
tual magnitude of this exponent becomes arbitrary and 
varies with the range of A9 we adopt for the power law 
fitting. This therefore casts doubts on the scaling argu- 
ments people conjure up to derive a and conclude that 



it plays the same universal role as the exponents in the 
critical phenomena [l7j . 

In order to study the dynamical behavior of Eq.(3), 
the actual dependence of £ on I is required. Note that 
our system size is growing in time. This is a moving 
boundary problem, and both the EW and KPZ equations 
fail to conserve the area. Instead of further complicating 
the existing models, we adopt the polar coordinates and 
modify the EW model to guarantee this conservation. 



OR 2 
~dt 



V ^§r + Z® [ r2 (t,o)-R 2 (t)} 



V a 



(3) 



where V a denotes the speed of fluid injection. Two phe- 
nomenological coefficients are introduced. First, the re- 
laxation coefficient v is the same as in the EW and KPZ 
models, which derives from the surface tension and fa- 
vors a smooth interface. The second term on the RHS is 
new. We ascribe it to the inertia effect which tends to 
magnify the existing crests by allowing them to advance 
faster in the next moment than the troughs. Therefore, 
we expect its coefficient £ to be proportional to the av- 
erage kinetic energy of the front, a (dR / dt) 2 ~ <rV a /t. 
It is interesting that Eq.(4), under the change of vari- 
able from 9 to x = R ■ 0, will generate a nonlinear term 
(dR/dx) 2 characteristic of the KPZ model. Dimensional 
analysis allows us to pinpoint the dependence of these 
two phenomcnological coefficients on the microscopic pa- 
rameters as: £ {f) = £ jt = KcrV a /r]r 2 t and v = Tr/rj 
where er, n, rj, r, T are the fluid density, permeability, vis- 
cosity, bead size, and dynamic surface tension. The rea- 
son why the contact angle did not come in is that we 
have limited the spacing between the two acrylic plates 
to be just a few beads to minimize its role. Note that the 
source term in Eq.(4) vanishes when we plug Eq.(l) into 
Eq.(4). Applying the orthonormal relation, we end up 
with a simple linear PDE for each of a; and bi . They can 
be easily solved to determine the evolution of the average 
noise spectrum Q as, 

Cf(t) = Cf(to)e- 2 ^ t - ta Ht/tor' (4) 
Substituting the new noise form into Eq.(3), we obtain 



^(^^^^^(M (5) 



N'sr 2 ,t^' 

where N' is the fluid-covered area divided by glass bead 
cross-section and Cz(irj) denotes the initial noise spectrum 
at in. Since there is no explicit random term in our PDE, 
this initial pattern cannot be taken as a perfect circle, 
otherwise the subsequent evolution will not be able to 
produce any roughness. But, we also need to make sure 
later that the final theory does not depend on the spe- 
cific choice of this initial time to and its corresponding 
pattern (i(to). We considered several possible situations 
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FIG. 3: Logarithmic plot of the time dependence of m defined 
in Eq.(9) for three invasion speeds, Vol = 2.50, V a 2 = 2 AO, 
and V a 3 = 2.19 cm 2 /s. The solid line in the inset indeed passes 
through the origin when extrapolated. 
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FIG. 4: Interface width versus time for invasion velocity 3.15 
cm 2 /s under three different A6. The forms of these three 
solid lines are determined by the model. The inset illustrate 
13 against A8 for invasion velocity 3.15 cm 2 /s. 



and generated the spectrum of the noise by analyzing the 
interface. The final form we adopt to describe the decay 
of noise with I is Q (to) — l/(s + I + s'l 2 ) where the value 
of s is found to fall between 0.1 and 1, and s' is around 
0.001. This corresponds to an initial shape in Fig. 2. 

Our experimental setup includes a thin cell (0.15 cm 
deep) formed by a horizontally positioned acrylic dish 
(radius=17 cm) and an acrylic top plate, a fluid injec- 
tion system, and a data recorder. The syringe injects the 
fluid vertically into the disorder media from the center 
of the dish with speeds ranging from 2.19 to 2.50 cm 2 /s. 
Glassbeads of diameter 500 /im are well tamped between 
the acrylic plates to make sure they will not move during 
the invasion. We use a high resolution (1 million pixels) 
CCD camera to monitor and record the interface, which 
is set vertically above the dish center. The camera mon- 
itors an area of 24 x 24 cm 2 plus the edge space, and the 
size of each pixel is about 450 /urn 2 . After the picture 
has been taken, a personal computer digitizes the video 
image and calculates the interface width. 

Experimentally the weighting of the noise spectrum is 
found to decrease with the Fourier angular momentum. 
This is evidenced by Fig. 1 , which shows the interface wig- 
gles progressively more severe with time. The high modes 
come in as more local structures sprout up at the newly 
gained surface. Lower modes, which correspond to more 
global structures, emerge earlier and thus have a longer 
history to evolve. According to the instability term in 
Eq. (4) , fluctuations are amplified via feedbacks on them- 
selves. Therefore, older structures tend to have larger 
amplitudes. The scale dependence of the width is shown 
in Fig. [5] The linear behavior of the width with a slope 
exactly equal to one in the asymptotic regime of small 
scales is consistent with we expected. Same features can 
be found in the bulk of non super-rough interface and 
many recent experiments 0, [l8l . Il9j ] . 

We use the parameters in Tab. II to estimate £' and v 
at about 1 and 0.036. Experimental values put them at 
about 0.8 ~ 1.1 (inset of Fig.4) and 10~ 2 ~ 10 _1 respec- 



TABLE II: The values of the relevant parameters in the ex- 
periment at room temperature. 



Parameter 


Value (unit) 


water-air surface tension (T) 


0.0727 N/m 


dynamic viscosity (77) 


1.002 x 10~ 3 Pa • s 


permeability (k) 


~ 10" 9 m 2 



tively. Other features are also verified experimentally, 
which include: (1) The instability coefficient £' is linearly 
proportional to the invasion velocity (Fig.4), (2) no de- 
pendence on V a for the relaxation coefficient v. Moreover, 
Fig. 5 finds A9 dependence for f3, as has been reported for 
a parallel injection^]. By use of Eq.(6), we can obtain 
the initial noise spectrum Q(to) by fitting the scale de- 
pendence of the interface width. Nicely, we deduce the 
same form at different times. Note that it is necessary to 
incorporate the possible symmetry of the initial configu- 
ration before being able to determine it directly from the 
spectrum in Fig. 2. We duly checked that the interface 
width was not sensitive to the exact choice of this ini- 
tial time to and shape of interface as long as they belong 
to the early stage when the effect of randomness is still 
small. As illustrated above, both a and j3 are smooth 
and varying functions of the scale. On the other hand, 
£' and v are good parameters to describe the scaling and 
dynamic behavior of our circular interface. 

Note that with our previous choice of a constant in- 
jection speed the advancing interface is decelerating. In 
order to check that our conclusions are not its artifacts, 
we repeated by accelerating the injection flux to attain a 
constant interface velocity. Now, £ = aKV^/rjr 2 becomes 
time- independent, V a is replaced by V^t, but parameter 
v is the same. The experimental results in Fig.6-(a),(b) 
turn out to be qualitatively the same as Fig. 3. We find 
that £ is still proportional to V' a with a slope that falls 
rightly within our estimated range an/rjr 2 = 10 3 ~ 10 4 
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FIG. 5: Interface width versus A9 at three different stages 
labelled by early time, transient, and saturation. The solid 
line comes from our formula. Invasion velocity is 5.1cm 2 /s. 
The inset demonstrates that the actual value of a varies with 
the range of AO we use to do the power law fitting. 
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(Fig.6-(c)) and v remains independent of V' a with the 
right magnitude (Fig.6-(d)). We believe that if the value 
of k can be more accurately measured, the correctness of 
an/rjr 2 shall be more solidly verified. We can now safely 
say that our dimensional analysis in the constant flow 
rate remains valid for the accelerated flow rate. 



In conclusion, we pointed out that the roughness of a 
circular interface is intrinsically different from the stan- 
dard planar injection. Major difference lies in the fact 
that the former interface grows in size. This results in a 
noise spectrum whose amplitudes decrease with momen- 
tum, in contrast to the white noise for a planar interface. 
One major consequence is the lack of power law in the 
scaling and dynamic behavior, which we believe should 
be common to all circular patterns [2uT | . Interestingly, we 
also find that the time exponent varies with the range of 
scale, as has been reported @ for a planar interface but 
not captured by existing models. All these features can 
be described by a linearized model we modified from the 
EW equation. Two phcnomenological coefficients are put 
forward to represent the competing trend of instability 
and smoothness. Their dependence on the microscopic 
parameters are determined from dimensional arguments, 
and estimated values conform with the experiments. 



FIG. 6: (a) Scale dependence of a constant advancing inter- 
face with Va = 4.01pixel 2 /s 2 . Solid line is the theoretical 
prediction. Because the flow rate is slow initially, more time 
is required for the designated early stage than in Fig. 4. In (b), 
a is shown to vary with the scale, (c) £ is plotted against V a '. 
The slope of the fitting line, 6840, falls within our predicted 
range, (d) Experimental values of v are independent of Va, 
consistent with our theory. 
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